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Abstract 

Context. The advent of deep multiwavelength extragalactic surveys has led to the necessity for advanced and fast 
methods for photometric analysis. In fact, codes which allow analyses of the same regions of the sky observed at 
different wavelengths and resolutions are becoming essential to thoroughly exploit current and future data. In this 
context, a key issue is the confusion (i.e. blending) of sources in low-resolution images. 

Ainns. We present t-phot, a publicly available software package developed within the astrodeep project, t-phot is 
aimed at extracting accurate photometry from low-resolution images, where the blending of sources can be a serious 
problem for the accurate and unbiased measurement of ffuxes and colours. 

Methods, t-phot can be considered as the next generation to tfit, providing significant improvements over and above 
it and other similar codes (e.g. CONVPHOt). t-PHOT gathers data from a high-resolution image of a region of the sky, 
and uses this information (source positions and morphologies) to obtain priors for the photometric analysis of the lower 
resolution image of the same field. t-PHOT can handle different types of datasets as input priors, namely i) a list of 
objects that will be used to obtain cutouts from the real high-resolution image; ii) a set of analytical models (as .fits 
stamps); hi) a list of unresolved, point-like sources, useful for example for far infrared wavelength domains. 

Results. By means of simulations and analysis of real datasets, we show that t-PHOT yields accurate estimations of 
fluxes within the intrinsic uncertainties of the method, when systematic errors are taken into account (which can be 
done thanks to a flagging code given in the output). t-PHOT is many times faster than similar codes like tfit and 
CONVPHOT (up to hundreds, depending on the problem and the method adopted), whilst at the same time being more 
robust and more versatile. This makes it an excellent choice for the analysis of large datasets. When used with the same 
parameter sets as for tfit it yields almost identical results (although in a much shorter time); in addition we show how 
the use of different settings and methods significantly enhances the performance. 

Conclusions, t-phot proves to be a state-of-the-art tool for multi wavelength optical to far-infrared image photometry. 
Given its versatility and robustness, t-phot can be considered the preferred choice for combined photometric analysis 
of current and forthcoming extragalactic imaging surveys. 

Key words. Galaxy, photometry, multiwavelength, software 


1. Introduction 


* T-PHOT is publicly available for downloading from 
www.astrodeep.eu/t-phot/ . 

** Scottish Universities Physics Alliance 


Combining observational data from the same regions of the 
sky in different wavelength domains has become common 
practice in the past few years (e.g. Agiieros et al. 2005 
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Figure 1. Schematic representation of the PSF-matched algorithm implemented in t-phot. Top: two objects are clearly 
detected and separated in the high-resolution detection image (blue line). The same two objects are blended in the low- 
resolution measurement image and have different colours (red line). Middle: the two objects are isolated in the detection 
image and are individually smoothed to the PSF of the measurement image, to obtain normalized model templates. 
Bottom: the intensity of each object is scaled to match the global profile of the measurement imag e. The scaling factors 
are found with a global minimization over the object areas. Image from De Santis et al. (2007). 


Obric et al.||2006l [Grogin et al.|[2QTTj and many oth- 
ers). However, the use of both space-based and ground- 
based imaging instruments, with different sensitivities, 
pixel scales, angular resolutions, and survey depths, raises 
a number of challenging difficulties in the data analysis pro¬ 
cess. 


In this context, it is of particular interest to obtain de¬ 
tailed photometric measurements for high-redshift galaxies 
in the near-infrared (NIR; corresponding to rest-frame op¬ 
tical) and far-infrared (FIR) domains. In particular, great 
attention must be paid to bandpasses containing spectral 
features which allow a thorough investigation of the sources, 
disentangling degenerate observational features, and ob¬ 
taining cruci al clues to the understanding of the gal actic 
physics (e.g. Daddi et al. 2004 Fontana et al. 2009). At 
2 : > 3, for example, photometry longward of if-band is 
needed to locate and measure the size of the Balmer break. 
A passive galaxy at z 6 (with the Balmer break ly¬ 
ing longward of the K-hand) can have i7-band and 3.6/im 
fluxes compatible, for example, with a star forming, dusty 
galaxy at 2 : 1, and TC—band photometry is necessary 

in order to disentangle the degeneracy. However, the lim¬ 
ited resolution of the ground based TC—band observations 
can impose severe limits on the reliability of traditional 
aperture or even Point Spread Function (PSF) fitting pho¬ 
tometry. In addition, IRAC photometry is of crucial im¬ 
portance so that reliable photometric redshifts of red and 
high- 2 : sources can be obtained, and robust stellar mass es¬ 
timates can be derived. 


To address this, a high-resolution image (HRI), for ex¬ 
ample obtained from the Hubble Space Telescope in the op¬ 
tical domain, can be used to retrieve detailed information 
on the positions and morphologies of the sources in a given 
region of the sky. Such information can be subsequently 
used to perform the photometric analysis of the lower reso¬ 
lution image (LRI), using the HRI data as priors. However, 
simply performing aperture photometry on the LRI at the 
positions measured in the HRI can be dramatically affected 


by neighbour contamination for reasonably sized apertures. 
On the other hand, performing source extraction on both 
images and matching the resulting catalogues is compro¬ 
mised by the inability to deblend neighbouring objects, 
and may introduce significant inaccuracies in the cross¬ 
correlation process. PSF-matching techniques that degrade 
high-resolution data to match the low-resolution data dis¬ 
card much of the valuable information obtained in the HRI, 
reducing all images to the “lowest common denominator” 
of angular resolution. Moreover, crowde d-field, PSF-fi tting 
photometry packages such as daophot ( Stetson|1987 ) per¬ 
form well if the sources in the LRI are unresolved, but are 
unsuitable for analysis of even marginally resolved images 
of extragalactic sources. 

A more viable approach consists of taking advantage of 
the morphological information given by the HRI, in order 
to obtain high-resolution cutouts or models of the sources. 
These priors can then be degraded to the resolution of 
the LRI using a suitable convolution kernel^ constructed 
by matching the PSFs of the HRI and of the LRI. Such 
low-resolution templates, normalized to unit flux, can then 
be placed at the positions given by the HRI detections, 
and the multiplicative factor that must be assigned to each 
model to match the measured flux in each pixel of the LRI 
will give the measured flux of that source. Such an ap¬ 
proach, although relying on some demanding assumptions 
as described in the following sections, has proven to be ef- 
flcient . It has been impl emented in such public codes as 
TFIT ( [Laidler et M^|2QQ7 ) and CONVPHOT ( |De Santis et al. 


previous 


2QQ7[ ), and h as already been utilized successfully in 
studies (e.g. Guo et ar]|2Q13 Galametz et al.||2Q13 ). 

In this paper we describe a new software package, T- 
PHOT, developed at INAF-OAR as part of the astrodeep 
projecfj^ The t-phot software can be considered a new, 
largely improved version of TFIT, supplemented with many 


^ astrodeep is a coordinated and comprehensive program of 
i) algorithm/software development and testing; ii) data reduc¬ 
tion/release, and iii) scientific data validation/analysis of the 
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of the features of CONVPHOT. Moreover, it adds many im¬ 
portant new options, including the possibility of adopting 
different types of priors (namely real images, analytical 
models, or point-sources). In particular, it is possible to use 
T-PHOT on FIR and sub-millimetric (sub-mm) datasets, as 
a competitive altern ative to the existing ded icated software 
such as FastPho t ( Bethermin et al.||2QlQ[ ) and DesPhot 
(Roseboom et al. ||2010| |Wang et al. j2014|). This makes T- 
PHOT a versatile tool, suitable for the photometric analysis 
of a very broad range of wavelengths from UV to sub-mm. 

T-PHOT is a robust and easy-to-handle code, with a 
precise structural architecture (a Python envelope calling 
C/C++ core codes) in which different routines are encap¬ 
sulated, implementing various numerical/conceptual meth¬ 
ods, to be chosen by simple switches in a parameter file. 
While a standard default “best choice” mode is provided 
and suggested, the user is allowed to select a preferred set¬ 
ting. 

One of the main advantages of t-phot is a significant 
saving of computational time with respect to both tfit and 
CONVPHOT (see Sect. |^. This has been achieved with the 
use of fast C modules and an efficient structural arrange¬ 
ment of the code. In addition to this, we demonstrate how 
different choices of parameters influence the performace, 
and can be optimized to significantly improve the final re¬ 
sults with respect to tfit, for example. 

The plan of the paper is as follows. Sect, [^provides a 
general introduction to the code, its mode of operation and 
its algorithms. In Section we discuss some assumptions, 
limitations and caveats of the method. Section presents 
a comprehensive set of tests, based on simulated and real 
datasets, to assess the performance of the code and to fully 
illustrate its capabilities and limitations. Section briefly 
discusses the computational performances of t-phot and 
provides some reference computational timescales. Finally, 
in Section the key features of t-phot are summarized, 
and outstanding issues and potential complications are 
briefly discussed. 


2. General description of the code 

As described above, t-phot uses spatial and morphologi¬ 
cal information gathered from a HRI to measure the fluxes 
in a LRI. To this end, a linear system is built and solved 
via matricial computing, minimizing the (in which the 
numerically determined fluxes for each detected source are 
compared to the measured fluxes in the LRI, summing the 
contributions of all pixels). Moreover, the code produces 
a number of diagnostic outputs and allows for an iterative 
re-calibration of the results. Figureshows a schematic de¬ 
piction of the basic PSF-matched fitting algorithm used in 
the code. 

As HRI priors t-phot can use i) real cutouts of sources 
from the HRI, ii) models of sources obtained with Galfit 
or similar codes, hi) a list of coordinates where PSF-shaped 
sources will be placed, or a combination of these three types 
of priors. 

For a detailed technical description of the mode of oper¬ 
ation of the code, we refer the reader to the Appendix and 
to the documentation included in the downloadable tarball. 
Here, we will briefly describe its main features. 


2.1. Pipeline 

The pipeline followed by t-phot is outlined in the flowchart 
given in Fig. The following paragraphs give a short de¬ 
scription of the pipeline. 


2.1.1. Input 

The input files needed by t-phot vary depending on the 
type(s) of priors used. 

If true high-resolution priors are used, e.g. for op- 
tical/NIR ground-based or IRAC measurements using 
Hubble Space Telescope (HST) cutouts, t-phot needs 

— the detection, high-resolution image (HRI) in .fits for¬ 
mat; 

— the catalogue of the sources in the HRI, obtained using 
SExtractor or similar codes (the required format is 
described in Appendix A); 

— the segmentation map of the HRI, in .fits format, 
again obtained using S Extractor or similar codes, 
having the value of the id of each source in the pix¬ 
els belonging to it, and zero everywhere else; 

— a convolution kernel LC, in the format of a .fits image 
or of a .txt file, matching the PSFs of the HRI and the 
LRI so that PSFlrj = K ^ PSFhri (* is the symbol 
for convolution). The kernel must have the HRI pixel 
scale. 


If analytical models are used as priors (e.g. Galfit 
models), T-PHOT needs 

— the stamps of the models (one per object, in .fits for¬ 
mat); 

— the catalogue of the models (the required format is de¬ 
scribed in Appendix A); 

— the convolution kernel K matching the PSFs of the HRI 
and the LRI, as in the previous case. 

If models have more than one component, one separate 
stamp per component and catalogues for each component 
are needed (e.g. one catalogue for bulges and one catalogue 
for disks). 

If unresolved, point-like priors are used, t-phot needs 

— the catalogue of positions (the required format is de¬ 
scribed in Appendix A); 

— the LRI PSF, in the LRI pixel scale. 


In this case, a potential limitation to the reliability of 
the method is given by the fact that the prior density usu¬ 
ally needs to be opt imized with resp ect to FIR/sub-mm 
maps, as discussed in| Shu etliL (|2Q15[ in preparation) and 
Elbaz et al. ( |2Q11 ) (see also |Wang et al.|2Q15||Bourne et al. 


2015, in preparation). The optimal nurnber of priors turns 


out to be around 50-75% of the numbers of beams in the 
map. The main problem is identifying which of the many 
potential priors from, for example, an HST catalogue one 
should use. This is a very complex issue and we do not 
discuss it in this paper. 

If mixed priors are used, t-phot obviously needs the 
input files corresponding to each of the different types of 
priors in use. 

Finally, in all cases t-phot needs 


deepest multi-wavelength cosmic surveys. For more information, 
visit http://astrodeep.eu . 


the measure LRI, background subtracted (see next para¬ 
graph), in .fits format, with the same orientation as 
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Figure 2. Schematic representation of the workflow in t-phot. 


the HRI (i.e. no rotation allowed); the pixel scale can be 
equal to, or an integer multiple of, the HRI pixel scale, 
and the origin of one pixel must coincide; it should be in 
surface brightness units (e.g. count s/s/pixel, or Jy/pixel 
for FIR images, and not PSF-filtered); 

— the LRI RMS map, in .fits format, with the same 
dimensions and WCS of the LRI. 

Table summarizes the input requirements for the dif¬ 
ferent choices of priors just described. 

All the input images must have the following key¬ 
words in their headers: CRPIXn, CRVALn, CDn_n, CTYPEn 
(n=l,2). 

2.1.2. Background subtraction 

As already mentioned, the LRI must be background sub¬ 
tracted before being fed to t-phot. This is of particular 
interest when dealing with FIR/sub-mm images, where the 
typical standard is to use zero-mean. To estimate the back¬ 
ground level in optical/NIR images, one simple possibility 
is to take advantage of the option to fit point-like sources 
to measure the flux for a list of positions chosen to fall 
within void regions. The issue is more problematic in such 
confusion-limited FIR images where there are no empty sky 
regions. In such cases, it is important to separate the fit¬ 
ted sources (those listed in the prior catalogue) from the 
background sources, which contribute to a flat background 
level behind the sources of interest. The priors should be 
chosen so that these two populations are uncorrelated. The 
average contribution of the faint background source pop¬ 
ulation can then be estimated for example by i) injecting 
fake sources into the map and measuring the average offset 
(output-input) flux, or ii) measuring the modal value in the 


residual image afte r a first pass through t-phot (see e.g. 


_ m _ 

Bourne et al.|2Q15[ in preparation). 


2.1.3. Stages 

T-PHOT goes through “stages”, each of which performs a 
well-defined task. The best results are obtained by per¬ 
forming two runs (“pass 1” and “pass 2”), the second us¬ 
ing locally registered kernels produced during the first. The 
possible stages are the following: 


— priors: creates/organizes stamps for sources as listed 
in the input priors catalogue(s); 

— convolve: convolves each high-resolution stamp with 
the convolution kernel K to obtain models (“tem¬ 
plates”) of the sources at LRI resolution. The templates 
are normalized to unit total flux. If the pixel scale of the 
images is different, transforms templates accordingly. 
Convolution is preferably performed in Fourier space, 
using fast FFTW3 libraries; however the user can choose 
to perform it in real pixel space, ensuring a more accu¬ 
rate result at the expense of a much slower computation. 

— positions: if an input catalogue of unresolved sources 
is given, creates the PSF-shaped templates listed in it, 
and merges it with the one produced in the convolve 
stage; 

— fit: performs the fitting procedure, solving the lin¬ 
ear system and obtaining the multiplicative factors to 
match each template flux with the measured one; 
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Table 1. The input files needed by t-phot for different settings. See text for details. 


— diags: selects the best fit!0 and produces the final for¬ 
matted output catalogues with fluxe s and errors, plus 
some other diagnostics, see Sect. |2.3| 

— dance: obtains local convolution kernels for the second 
pass; it can be skipped if the user is only interested in 
a single-pass run; 

— plotdance: plots diagnostics for the dance stage; it can 
be skipped for any purpose other than diagnostics; 

— archive: archives all results in a subdirectory whose 
name is based on the LRI and the chosen fitting method 
(to be used only at the end of the second pass). 

The exact pipeline followed by the code is specified by 
a keyword in the input parameter file. See also Appendix 
A for a more detailed description of the whole procedure. 


2.1.4. Solution of the linear system 


The search for the LRI fluxes of the objects detected in the 
HRI is performed by creating a linear system 

N 

= Y^J2FiPi{m,n) ( 1 ) 

m,n m,n i 

where m and n are the pixel indexes, / contains the pixel 
values of the fluxes in the LRI, Pi is the normalized flux of 
the template for the i-th objects in the (region of the) LRI 
being fitted, and Fi is the multiplicative scaling factor for 
each object. In physical terms, Fi represents the flux of each 
object in the LRI (i.e. it is the unknown to be determined). 

Once the normalized templates for each object in the 
LRI (or region of interest within the LRI) have been gener¬ 
ated during the convolve stage, the best fit to their fluxes 
can be simultaneously derived by minimizing a statistic, 

2 \Em,nIi'm,n)-M{m,n)f 

X = - ^-7-7- (2) 

a{m,n) 

where m and n are the pixel indexes, 

N 

M{m,n) ='^FiPi{m,n) (3) 

i 


and a is the value of the RMS map at the (m, n) pixel 
position. 

The output quantities are the best-fit solutions of the 
minimization procedure, i.e. the Fi parameters and their 
relative errors. They can be obtained by resolving the linear 
system 



(4) 


^ Each source is fitted more than once if an arbitrary grid is 
used, as in the standard tfit approach. 


for i = 0,1,..., N. 

In practice, the linear system can be rearranged into a 
matrix equation, 

AF = B (5) 

where the matrix A contains the coefficients PiP^ja^^ F 
is a vector containing the fluxes to be determined, and B 
is a vector given by liPija^ terms. The matrix equation is 
solved via one of three possible methods as described in the 
next subsection. 


2.1.5. Fitting options 

T-PHOT allows for some different options to perform the fit: 


three different methods for solving the linear system 
are implemented, namely, the LU method (used by de¬ 
fault in tfit), the Cholesky method, and the Iterative 
Biconjugate Gradient method (used by default in CON- 
VPHOT; for a r eview on me thods to solve sparse linear 
systems see e.g. |Davis|2QQ6 ). They yield similar results, 
although the LU method is slightly more stable and 
faster; 

a threshold can be imposed so that only pixels with a 
flux higher than th is leve l will be used in the fitting 
procedure (see Sect. 4.1.4); 

sources fitted with a large, unphysical negative flux 
{fmeas < —3(7, where (7 is their nominal error, see be¬ 
low) can be excluded from the fit, and in this case a new 
fitting loop will be performed without considering these 
sources. 


The fit can be performed i) on the entire LRI as a whole, 
producing a single matrix containing all the sources (this is 
the method adopted in CONVPHOt); ii) subdividing the LRI 
into an arbitrary grid of (overlapping) small cells, perfom- 
ing the fit in each of such cells separately, and then choosing 
the best fit for each source, using some convenient criteria 
to select it (because sources will be fitted more than once if 
the cells overlap; this is the method adopted in tfit); iii) 
ordering objects by decreasing flux, building a cell around 
each source including all its potential contaminants, solv¬ 
ing the problem in that cell and assigning to the source the 
obtained flux {cells-on-objects method; see the Appendix 
for more details). 

While the first method is the safest and more accurate 
because it does not introduce any bias or arbitrary modifi¬ 
cations, it may often be unfeasible to process at once large 
or very crowded images. Potentially large computational 
time saving is possible using the cells-on-objects method, 
depending on the level of blending/confusion in the LRI; if 
it is very high, most sources will overlap and the cells will 
end up being very large. This ultimately results in repeating 
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Figures. Example of the results of a standard t-phot run using extended priors. Left to right: HRI (FWHM=0.2”), 
LRI (FWHM=1.66”), and residuals image for a simulated dataset. LRI and residual image are on the same greyscale. 



Figure 4. Example 
and residuals image 


of the results of a standard t-phot run using point-sou rce pr iors. Left to right: LRI (FWHM=25”) 


(same greyscale) for a simulated dataset. See also Sect. 4.1.2 



Figure 5. Example of the results of a standard t-phot run using analytical priors. Left to right: CANDELS COSMOS 
iL-band (HRI), R-band (LRI) and residuals image obtained using Galfit two-component models. LRI and residual 
image are on the same greyscale. 


many times the fit on regions with dimensions comparable 
to the whole image (a check is implemented in the code, 
to automatically change the method from cells-on-objects 
to single fit if this is the case). If the confusion is not dra¬ 
matic, a saving in computational time up to two orders of 
magnitude can be achieved. The results obtained using the 


cells-on-objects method prove to be virtually identical to 
those obtai ned with a single fit on the whole image (see 
Sect. 4.1.2). On the other hand, using the arbitrary cells 
method is normally the fastest option, but can introduce 
potentially large errors to the flux estimates owing to wrong 
assignments of peripheral fiux from sources located outside 
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a given cell to sources within the cell (again, see Sect. 4.1.2 
and the Appendix B). 


-h4 if the source is at the border of the image (i.e. its seg¬ 
mentation reaches the limits of the HRI pixels range). 


2.1.6. Post-fitting stages: kernel registration 

After the fitting procedure is completed, t-phot will pro¬ 
duce the f inal output catalogues and diagnostic images (see 
Sect. |2.3[ ). Among these, a model image is obtained by 
adding all the templates, scaled to their correct total fiux 
after fitting, in the positions of the sources. This image will 
subsequently be used if a second pass is planned; during a 
stage named dance, a list of positional shifts is computed, 
and a set of shifted kernels are generated and stored. The 
dance stage consists of three conceptual steps: 

— the LRI is divided into cells of a given size (specified 
by the keyword dzonesize) and a linear Ax, shift is 
computed within each cell, cross-correlating the model 
image and the LRI in the considered regiorjj 

— interpolated shifts are computed for the regions where 
the previous registration process gives spuriously large 
shifts, i.e. above the given input threshold parameter 
maxshift; 

— the new set of kernels is created using the computed 
shifts to linearly interpolate their positions, while cata¬ 
logues reporting the shifts and the paths to kernels are 
produced. 

2.1.7. Second pass 

The registered kernels can subsequently be used in the sec¬ 
ond pass run to obtain more astrometrically precise results. 
T-PHOT automatically deals with them provided the correct 
keyword is given in the parameter file. If unresolved priors 
are used, the list of shifts generated in the dance stage will 
be used by the positions routine during the second pass 
to produce correctly shifted PSFs and generate new tem¬ 
plates. 

2.2. Error budget 

During the fitting stage, the covariance matrix is con¬ 
structed. Errors for each source are assigned as the square 
root of the diagonal element of the covariance matrix rela¬ 
tive to that source. It must be pointed out that using any 
cell method for the fitting rather than the single fitting 
option will affect this uncertainty budget, since a different 
matrix will be constructed and resolved in each cell. 

It is important to stress that this covariance error bud¬ 
get is a statistical uncertainty, relative to the RMS fluc¬ 
tuations in the measurement image, and is not related to 
any possible systematic error. The latter can instead be 
estimated by flagging potentially problematic sources, to 
be identified separately from the fitting procedure. There 
can be different possible causes for systematic offsets of the 
measured flux with respect to the true flux of a source. 
T-PHOT assigns the following flags: 

-hi if the prior has saturated or negative flux; 

— +2 if the prior is blended (the check is performed on the 
segmentation map); 

3 pp^ direct cross-correlations are implemented, the lat¬ 
ter being the preferred default choice because it gives more pre¬ 
cise results at the expense of a slightly slower computation. 


2.3. Description of the output 

T-PHOT output files are designed to be very similar in for¬ 
mat to those produced by tfit. They provide 

— a “best” catalogue containing the following data, listed 
for each detected source (as reported in the catalogue 
file header): 

— id; 

— X and y positions (in LRI pixel scale and reference 
frame, FITS convention where the first pixel is cen¬ 
tred at 1,1); 

— id of the cell in which the best fit has been obtained 
(only relevant for the arbitrary grid fitting method); 

— X and y positions of the object in the cell and dis¬ 
tance from the centre (always equal to 0 if the cells- 
on-objects method is adopted); 

— fitted fiux and its uncertainty (square root of the 
variance, from the covariance matrix). These are the 
most important output quantities; 

— fiux of the object as given in the input HRI catalogue 
or, in the case of point-source priors, measured flux 
of the pixel at the x, y position of the source in the 
LRI; 

— fiux of the object as determined in the cutout stage 
(it can be different to the previous one, e.g. if the seg¬ 
mentation was dilated); in the case of point-sources 
priors, measured fiux of the pixel at the x, y position 
of the source in the LRI; 

— flag indicating a possible bad source as described in 
the previous subsection; 

— number of fits for the object (only relevant for arbi¬ 
trary grid fitting method, 1 in all other cases). 

— id of the object having the largest covariance with 
the present source; 

— covariance index, i.e. the ratio of the maximum co- 
variance to the variance of the object itself; this num¬ 
ber can be considered an indicator of the reliability 
of the fit, since large covariances often indicate a 
possible systematic offset in t he me asured flux of 
the covarying objects (see Sect. |4.1.^ ). 

— two catalogues reporting statistics for the fitting cells 
and the covariance matrices (they are described in the 
documentation); 

— the model .fits image, obtained as a collage of the 
templates, as already described; 

— a diagnostic residual .fits image, obtained by sub¬ 
tracting the model image from the LRI; 

— a subdirectory containing all the low-resolution model 
templates; 

— a subdirectory containing the covariance matrices in 
graphic (.fits) format; 

— a few ancillary files relating to the shifts of the kernel 
for the second pass and a subdirectory containing the 
shifted kernels. 

All fluxes and errors are output in units consistent with 
the input images. 

Figuresandshow three examples of t-phot appli¬ 
cations on simulated and real data, using the three different 
options for priors. 
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3. Assumptions and limitations 

The PSF-matching algorithms implemented in t-phot and 
described in the previous section are prone to some assump¬ 
tions and limitations. In particular, the following issues 
must be pointed out. 

i) The accuracy of the results strongly depends on the 
reliability of the determined PSFs (and consequently of the 
convolution kernel). An error of a few percentage points in 
the central slope of the PSF light profile might lead to non- 
negligible systematical deviations in the measured fluxes. 
However, since the fitting algorithm minimizes the residu¬ 
als on the basis of a summation over pixels, an incorrect 
PSF profile will lead to characteristic positive and nega¬ 
tive ring-shaped patterns in the residuals (see Fig. [^, and 
to some extent the summation over pixels will compensate 
the global fiux determination. 



Figure 6. Typical patterns in a residual image created by 
T-PHOT, caused by inaccurate PSF/kernel determination. 
In this case, the ring-shaped shadow surrounding a bright 
central spot is due to an underestimation of the central 
peak of the LRI PSF, which causes an overestimation of 
the fit in the outskirts while leaving too much light in the 
centre. 


ii) When dealing with extended priors, it is assumed 
that the instrinsic morphology of the objects does not 
change with the wavelength. Of course, this is usually not 
the case. The issue is less of a problem when dealing with 
FIR images, in which the morphological features of the pri¬ 
ors are unresolved by the low-resolution PSF. On the other 
hand, in the optical and NIR domains this problem may 
be solved by the use of multicomponent analytical models 
as priors. In this approach, each component should be fit¬ 
ted independently, thus allowing the ratio between bulge 
and disk components to vary between the HRI and LRI. A 
clear drawback of this approach is that any failure of the fit 
due to irregular or difficult morphological features (spiral 
arms, blobs, asymmetries, etc.) would be propagated into 
the LRI solution. This functionality is already implemented 
in T-PHOT and detailed testing is ongoing. 

hi) As explained in Sect. 2.2, t-phot flags priors that 
are likely to be flawed: sources too close to the borders 
of the image, saturated objects, and most notably blended 
priors. The assumption that all priors are well separated 
from one another is crucial, and the method fails when this 
requirement is not accomplished. Again, this is crucial only 
when dealing with real priors, while analytical models and 
unresolved priors are not affected by this limitation. 


iv) As anticipated in Sect. 2.1.1, FIR images can suffer 
from an “overfitting” problem, due to the presence of too 
many priors in each LRI beam if the HRI is deeper than the 
LRI. In this case, a selection of the priors based on some ad¬ 
ditional criteria (e.g. fiux predition from SED fitting) might 


be necessary to avoid catastrop hic outcomes (see also [Wang 


et al. 1(2^5 |Bourne et al.||2Q15 in preparation). 


4. Validation 


To assess the performance of t-phot we set up an extensive 
set of simulations, aimed at various different and comple¬ 
mentary goals. 


We used SkyMaker ( |Bertin||2QQ9[ ), a public software 
tool, to build synthetic .fits images. The code ensures di¬ 
rect control on all the observational parameters (the mag¬ 
nitude and positions of the objects, their morphology, the 
zero point magnitude, the noise level, and the PSF). Model 
galaxies were built by summing a de Vaucouleurs and an 
exponential light profile in order to best mimic a realis¬ 
tic distribution of galaxy morphologies. These models were 
generated using a variety of bulge-to-total light ratios, com¬ 
ponent sizes, and projection angles. 

All tests were run using ideal (i.e. synthetic and sym¬ 
metric) PSFs and kernels. 

Moreover, we also performed tests on real datasets taken 
from the CANDELS survey (in these cases using real PSFs). 

Some of the tests were performed using both t-phot 
and TFIT, to cross-check the results, ensuring the perfect 
equivalency of the results given by the two codes when used 
with the same parameter sets, and showing how appropriate 
settings of the t-phot parameters can ensure remarkable 
improvements. 

For simplicity, here we only show the results from a re¬ 
stricted selection of the test dataset, which are representa¬ 
tive of the performance of t-phot in standard situations. 
The results of the other simulations resemble overall the 
ones we present, and are omitted for the sake of concise- 


4.1. Code performance and reliability on simulated images 

4.1.1. Basic tests 

As a first test, we checked the performance of the basic 
method by measuring the flux of two PSF-shaped synthetic 
sources, with varying separation and signal-to-noise ratios. 
One hundred realizations with different noise maps of each 
parameter set were prepared, and the averages on the mea¬ 
sured fluxes were computed. The aims of this test were 
twofold: on the one hand, to check the precision to which 
the fitting method can retrieve true fluxes in the simplest 
possible case - two sources with ideal PSF shape; on the 
other hand, to check the reliability of the nominal error 
budget given by the covariance matrix, comparing it to the 
real RMS of the 100 measurements. Figure shows three 
examples of the set-up and the results of this test. Clearly, 
in both aspects the results are reassuring: the average of the 
100 measurements (red diamonds) is always in very good 
agreement with the true value, with offset in relative er¬ 
ror always well under the 1/{S/N)centr limit {{S/N)centr is 
the value of the signal-to-noise ratio in the central pixel of 
the source, corresponding to roughly one third of the total 
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Figure 7. Accuracy check on idealized PSF-shaped objects. 100 realizations of the same image containing two PSF- 
shaped objects at varying positions and signal-to-noise ratios have been produced and the fluxes have been measured 
with T-PHOT. In each row, the left image shows one of the 100 realizations with the largest considered separation (10 
pixels). On the right, the first panel refers to the central object, and the second (on the right) to the shifted object; 
central signal-to-noise {S/Ncentr) ratios are, from top to bottom, 100, 100, 100 for the first source and 100, 10, 1 for the 
second source. In each panel, as a function of the separation interval between the two sources, the faint grey points show 
each of the 100 flux measurements (in relative difference with respect to the true input flux), the red diamonds are the 
averages of the 100 measurements, the red crosses show the nominal error given by the covariance matrix in t-phot, 
and the green dots the standard deviation of the 100 measurements. See text for more details. 


S/N); and the nominal error (red crosses) given by the co- 
variance matrix is always in good agreement with the RMS 
of the 100 measurements (red circles). 

When dealing with extended objects rather than with 
point-like sources, one must consider the additional prob¬ 
lem that the entire profile of the source cannot be measured 


exactly because the segmentation is limited by the lowest 
signal-to-noise isophote. The extension of the segmentation 
therefore plays a crucial role and defining it correctly is 
a very subtle problem. Simply taking the isophotal area 
as reported by SExtractor as ISOAREA often underesti¬ 
mates the real extension of the objects. Accordingly, the 
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Figure 8. Effects of different segmentation areas on the measured flux of two isolated objects with identical flux and 
signal-to-noise ratio, at two possible separations of 40 and 120 pixels. Each panel shows the flux error in one of the objects 
at each separation distance. The shades and dimensions of the dots is a function of the radius of the segmentation, with 
darker and smaller dots corresponding to smaller segmentations. See text for more details. 


segmentation of the sources should somehow be enlarged 
to include the faint wings of sources. To this aim, specific 
software called Dilate was developed at OAR and used 
in the CANDELS pipeline for the photometric analysis o f 


GOODS-S and UDS IRAC data jGalametz et al. 


€ 


2013). 


Dilate enlarges the segmentation by a given factor, de- 
pending on the original area; it has proven to be reasonably 
robust in minimizing the effects of underestimated segmen¬ 
tated areas. 

Figure shows the effects of artificially varying the di¬ 
mensions of the segmentation relative to two bright, ex¬ 
tended and isolated sources in a simulated HRI, on the flux 
measured for that source in a companion simulated LRI. It 
is important to note how enlarging the segmented area nor¬ 
mally results in larger measured fluxes, because more and 
more light from the faint wings of the source are included in 
the fit. However, beyond a certain limit the measurements 
begin to lose accuracy owing to the inclusion of noisy, too 
low signal-to-noise regions (which may cause a lower flux 
measurement). 

In principle, using extended analytical models rather 
than real high-resolution cutouts should cure this problem 
more efficiently, because models have extended wings that 
are not signal-to-noise limited. Tests are ongoing to check 
the performance of this approach, and will be presented in 
a forthcoming paper. 


4.1.2. Tests on realistic simulations 


The next tests were aimed at investigating less idealized sit¬ 
uations, and have been designed to provide a robust analy¬ 
sis of the performance of the code on realistic datasets. We 
used the code Gen Cat (Schreiber et al.||2Q 15, 

it he 


__ m prepara¬ 
tion) to produce mock catalogues of synthetic extragalactic 


sources, with reasonable morphological features and flux 
distributior]^ Then, a set of images were produced using 


^ GenCat is another software package developed within the 
ASTRODEEP project. It uses GOODS-S CANDELS statistics to 
generate a realistic distribution of masses at all redshifts, for 


such catalogues as an input for SkyMaker. A “detection” 
HRI mimicking an HST H band observation (FWHM = 
0.2") was generated from the GenCat catalogue using out¬ 
put parameters to characterize the objects’ extended prop¬ 
erties. Then a set of measure LRIs were produced: the first 
was populated with PSF-shaped sources, having FWHM = 
1.66” (the typical IRAC-chl/ch2 resolution, a key applica¬ 
tion for t-phot), while other LRIs were created from the 
input catalogue, mimicking different ground-based or IRAC 
FWHMs. Finally, we created another HRI catalogue remov¬ 
ing all of the overlapping source^ This “non-overlapping” 
catalogue was used to create parallel detection and mea¬ 
surement images in order to obtain insight into the compli¬ 
cations given by the presence of overlapping priors. In all 
these images, the limiting magnitude was set equal to the 
assigned zero point, so that the limiting flux at la is 1. In 
addition, the fits were always performed on the LRI as a 
whole, if not otherwise specified. 

Figure shows the results relative to the first test, 
i.e. the fit on the image containing non-overlapping, PSF- 
shaped sources, with a “perfect” detection (i.e. the priors 
catalogue contains all sources above the detection limit), 
obtained with a single fit on the whole image. The figure 
shows the relative error in the measured flux of the sources, 
ifmeas “ ftrue) ! ftrue, vcrsus the log of the real input flux 


two populations of galaxies (active and passive), consistently 
with observed mass functions. All the other physical proper¬ 
ties of the mock galaxies are then estimated using analytical 
recipes from literature: each source is assigned a morphology 
(bulge-to-total ratio, disk and bulge scale lengths, inclination 
etc.), star formation rate, attenuation, optical and infrared rest- 
frame, and observed magnitudes. Each source is finally assigned 
a sky-projected position mimicking the clustering properties of 
the real CANDELS data. 

^ We proceeded as follows. First, we created a “true” segmen¬ 
tation image using the input catalogue and assigning to each 
object all the pixels in which the flux was 1.005 x fbackground ■ 
Then, starting from the beginning of the list, we included each 
source in the new catalogue if its segmented area did not overlap 
the segmented area of another already inserted source. 
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Figure 9. Accuracy of the flux determination in a simulation containing non-overlapping, PSF-shaped sources and 
“perfect” detection. Relative measured flux difference {fmeas — ftrue)/ftrue is plotted versus logarithm of the input 
flux ftrue', for ^ simulated image populated with PSF-shaped sources (FWHM=1.66”). Each dot corresponds to a single 
source, with different symbols and colours referring to various diagnostics as explained in the legend and in the colourbar. 
The black solid line is the average in bins, the yellow shade is the standard deviation. The vertical dashed line shows 
the limiting flux at la, / = 1. The inner panel shows a magnification of the brighter end of the distribution. The fit was 
performed on the whole image at once. See text for more details. 


ftrue] fhe different symbols refer to the flag assigned to each 
object, while the colour is a proxy for the covariance index. 

In this case, the only source of uncertainty in the mea¬ 
surement is given by the noise fluctuations, which clearly 
become dominant at the faint end of the distribution. 
Looking at the error bars of the sources, which are given 
by the nominal error assigned by T-PHOT from the covari¬ 
ance matrix, one can see that almost all sources have mea¬ 
sured flux within 2a from their true flux, with only strongly 
covariant sources (covariance index 1, greener colours) 
having \ fmeas - ftruel/ftrue > IcT- The Only noticeable ex¬ 
ceptions are sources that have been flagged as potentially 
unreliable, as described in Sect. |2.2[ We also note how the 
average A/// (solid black line) is consistent with zero down 
to ftrue = SjN ^ 0.63. 

Figure pT| shows the analysis of a case study in which the 
fluxes of a clump of highly covariant objects are measured 
with poor accuracy, and some of the nominal uncertainties 
are underestimated: a very bright source (ID 3386, m^rue = 
21.17) shows a relative difference {fmeas-ftrue)/ftrue > 3a. 
To cast light on the reason for such a discrepancy, the region 
surrounding the object was replicated 100 times with dif¬ 
ferent noise realizations, and the results were analysed and 
compared. The upper panels show (left) one of the 100 mea¬ 
surement images and (right) the position of all the sources 
in the region (many of which are close to the detection 
limit). The colour code refers to the covariance index of the 
sources. The bottom left panel shows the relative error in 


the measured flux for all the sources in the region, with the 
inner panels showing magnifications relative to the object 
ID 3386 and to the bunch of objects with rn^me ^ 26.5. 
Looking at the colours of their symbols, many objects in 
the region turn out to be strongly covariant. Indeed, while 
the bluer sources in the upper part of the region all have 
covariant indexes lower than 0.5, the greener ones in the 
crowded lower part all have covariance index larger than 
I (indeed larger than 2 in many cases). This means that 
their flux measurements are subject to uncertainties not 
only from noise fluctuations, but also from systematic er¬ 
rors due to their extremely close and bright neighbours. As 
clearly demonstrated here, the covariance index can give a 
clue to which measurements can be safely trusted. 

The bottom right panel gives the sum of the measured 
fluxes of all sources in each of the 100 realizations (the blue 
line is the true total flux and the red line is the mean of 
the 100 measured total fluxes). It can be seen that the total 
flux measured in the region is always consistent with the 
expected true one to within ^ 1% of its value. 

Although it is not possible to postulate a one-to-one re¬ 
lation (because in most cases sources having a large covari¬ 
ance index have a relatively good flux estimate, see Fig.[^, 
the bottom line of this analysis is that the covariance index, 
together with the flagging code outputted by t-phot, can 
give clues about the reliability of measured flux, and should 
be taken into consideration during the analysis of the data. 
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Figure 10. Analysis of a small region including a strongly covarying group of sources. Upper left panel: one of the 100 
realizations with different noise maps of the region. Upper right panel: true spatial position of all the sources in the 
region (the colour of the dots refer to the covariance index of the sources, as indicated in the colourbar, while their size 
is proportional to their true flux). Bottom left panel: relative deviation of measured flux from the true flux for each 
source in the region, as a function of their true magnitude (big dots show the average relative deviation, and their colours 
refer to their covariance index as in the previous panel; green squares show the nominal uncertainty given by t-phot, 
to be compared with the RMS of the distribution of the 100 measurements (diamonds); small grey dots are the 100 
measurements. The insets show magnifications of regions of interest). Bottom right panel: each dot shows the sum of the 
measured fluxes for each of the 100 realizations, and the average of this sum (red line) to be compared with the true sum 
(blue line), showing that an overall consistency is guaranteed by the method. See text for more details. 


Measurements relative to sources having a covariance index 
larger than 1 should be treated with caution. 

In a subsequent more realistic test, we considered ex¬ 
tended objects (including morphologies of objects from 
the GenCat catalogue, using FWHMi:/i?/=0.2” and 
FWHMLi 7 / = 1.66” and imposing VCltrue,LRI=^true,HRI = 
GenCat for simplicity) and allowed for overlapping 
priors. To be consistent with the standard procedure 
adopted for real images, for this case we proceeded by pro¬ 
ducing an S Extractor catalogue and segmentation map, 
which were then spatially cross-correlated with the “true” 
input catalogue. The results for this test are shown in Fig. 
[TH Even in this much more complex situation, the results 
are reassuring; there is an overall good agreement between 
measured and input fluxes for bright (log S/N > 1) sources, 
with only a few flagged objects clearly showing large devi¬ 
ations from the expected value, and a reasonably good av¬ 
erage agreement down to log S'/A" = 0. However, all bright 
fluxes are measured 5% fainter than the true values (see 
the inner box in the same figure); this is very likely the effect 
of the limited segmentation extension, as already discussed 
in the previous section. On the other hand, faint sources 
tend to have systematically overestimated fluxes, arguably 


because of contamination from undetected sources. To con¬ 
firm this, we focus our attention on a single case study (the 
source marked as ID 720) which shows a large discrepancy 
from its true flux, but has a relatively small covariance in¬ 
dex. An analysis of the real segmentation map shows how 
the detected object is actually a superposition of two dif¬ 
ferent sources that have been detected as a single one, so 
that the measured flux is of course higher than expected. 
One should also note that the uncertainties on the measured 
fluxes are smaller in this test, because there are fewer priors 
(only the ones detected by SExtractor are now present), 
implying a lower rank of the covariance matrix and a lower 
number of detected neighbours blending in the LRI. This 
causes a global underestimation of the errors. 

To check the performance of t-phot at FIR wave¬ 
lengths, we also run a test on a simulated Herschel SPIRE 
250 yum image (FWHM=25”, 3.6” pixel scale). The simu¬ 
lated image (shown alongside with the obtained residuals 
in Fig. mimics real images from the GOODS-Herschel 
program, the deepest Herschel images ever obtained. This 


image was p roduced with the technique presented in Leiton 


et al. (2015); we first derived (predicted) flux densities for 
all the 24 /im detections (7^24^^ > 20/iJy) in GOODS- 
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Figure 11. Accuracy of the flux determination in a simulation containing extended objects, overlapping priors, and 
SExtractor detection. Top: relative flux difference {fmeas — ftrue)/ftrue versus logarithm of the input flux ftme for 
a simulated image populated with extended sources (FWHM=1.66”). Symbols and colours are as in Fig. The inner 
panel shows a magnification of the brighter end of the distribution. The outlier marked with the open black circle, 
ID=720, is shown in the bottom panel: left to right, HRI (FWHM=0.2”), LRI, SExtractor segmentation map and 
“true” segmentation map. The green circles show the object detected via SExtractor, while the blue cross shows its 
“true” position. See text for more details. 


North, which are dependent on their redshift and flux den¬ 
sities at shorter wavelengths, and then we injected these 
sources into the real noise maps from GOODS-Herschel 
imaging. Additional positional uncertainties, typical ly 0.5^^ 
were a l so app lied to mimic real images. As shown in |Leiton| 


et al.| ( |2Q15| ), these simulated images h ave similar pixel 


value distributions to real images (see also Wang et al. 2015 
for more details). For this test, t-phot was run using the 
list of all the 24 /im sources as unresolved priors. The results 
of the test are plotted in Fig.[^ and they show that even in 
this case t-phot can recover the input fluxes of the sources 
with great statistical accuracy (the mean of the relative de¬ 


viation from the expected measurements, i.e. the black solid 
line in the plot, is consistent with zero down to the faintest 
fluxes). The results are equivalent to those obtained on the 
same datasets with other public software specifi cally devel 


oped for FI R photometry, such as FastPhot (Bethermin 




4.1.3. Testing different fitting options: cell dimensions 

We then proceeded to check the performance of the dif¬ 
ferent fitting techniques that can be used in T-PHOT. To 
this aim, we repeated the test on the 1.66” LRI with ex- 
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Figure 12. Accuracy of the flux determination in a FIR-like simulation (Herschel SPIRE 250 yum, FWHM=25”, 3.6” 
pixel scale), using unresolved priors. The symbols have the same meaning as in Fig. See text for more details. 
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Figure 13. Accuracy of the flux determination. Top panel: for the same simulation described in Fig. 11, the histograms 


show the results for three different fitting methods: regular grid 100x100 pixels (standard tfit approach), regular grid 
200 x200, single fit on the whole image. The small boxes show the extended wings of the histograms, magnified for better 
viewing. The accuracy increases by enlarging the cells and, reaches the best result with the single fit on the whole image. 
Bottom panel: the histogram shows the relative measured flux difference between the single fit on the whole image and 
the cells-on-objects method. Differences above 1% are very rare. 


tended priors and SExtractor priors, described in Sect. 
|4.1.2[ with different fitting methods: using a regular grid of 
cells of 100 X 100 pixels, a regular grid of cells of 200 x 200 
pixels, and the cells-on-objects method, comparing the re¬ 
sults with those from the fit of the whole image at once. 
The results of the tests are shown in Figs. [T^ and [T^ The 


first figure compares the distributions of the relative errors 
in measured flux for the runs performed on the 100 x 100 
pixels grid, on the 200 x 200 pixels grid, and on the whole 


image at once. Clearly, using any regular grid of cells wors¬ 
ens the results, as anticipated in Sect. |2.1.'5l Enlarging the 
sizes of the cells improves the situation, but does not com¬ 
pletely solve the problem. We note that the adoption of an 
arbitary grid of cells of any dimension in principle is prone 
to the introduction of potentially large errors, because (pos¬ 
sibly bright) contaminating objects may contribute to the 
brightness measured in the cell, without being included as 
contributing sources. A mathematical sketch of this issue is 
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Figure 14. Accuracy of the flux determination. For the same simulation described in Fig.[^ the plots show the results 
for four different fitting methods. Top panel, left to right: HRI (FWHM=0.2”), LRI (FWHM=1.66”), residuals using a 
regular grid of 100 x 100 pixels cells (standard tfit approach), a regular grid of 200 x 200 pixels cells, a single fit on the 
whole image, and the cells-on-objects method. The spurious fluctuations in the last two panels are due to segmentation 
inaccuracies, as in Fig.[^ Bottom panels, left to right and top to bottom: relative measured flux differences with respect 
to true fluxes, same order as above. The values of the covariance index are different in each case because of the varying 
sizes of the cells (and therefore of the relative matrix). 


explained in the Appendix B (see also Sect. 4.2). The second 
histogram compares the differences between the fit on the 
whole image and the fit with the cells-on-objects method. 
Almost all the sources yield identical results with the two 
methods, within {fmeas-ftrue)/ ftrue < 0.001, which proves 
that the cells-on-objects method can be considered a reli¬ 
able alternative to the single-fit method. Finally, Figure 
compares the HRI, the LRI, and the residual images ob¬ 
tained with the four methods and their distributions of rel¬ 
ative errors, showing quantitatively the difference between 
the analyzed cases. 


and the equivalence of the results obtained with the single- 
fit and the cells-on-objects techniques, are evident. As al¬ 
ready anticipated, one should bear in mind that the cells- 
on-objects method is only convenient if the overlapping of 
sources is not dramatic, as in ground-based optical obser¬ 
vations. For IRAC and FIR images, on the other hand, 
the extreme blending of sources would cause the cells to be 
extended over regions approaching the size of the whole im¬ 
age, so that a single fit would be more convenient, although 
often still CPU-time consuming. 


In summary, it is clear that an incautious choice of cell 
size may lead to unsatisfactory and catastrophic outcomes. 
On the other hand, the advantages of using a single fit. 
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Threshold 0.0 


Threshold 0.5 


Threshold 0.9 



Figure 15. Effects of threshold fitting (Sect. 4.1.4). Mean relative error (black line) and standard deviation (yellow 
shaded area) for three simulations with different threshold values (0.0, 0.5, 0.9). Only pixels with normalized fiux higher 
than the threshold values are included in the fit. Larger threshold values result in more accurate measurements for faint 
sources, at the expense of a systematic underestimation of the fiux for brighter ones. 



X SEx MAG BEST 
o SEx MAGAPER 
T-PHOT 



''true 


Figure 16. Top: measured magnitude differences {rrimeas — '^true) versus “true” input magnitudes rutrue^ for two simu¬ 
lated images populated with extended sources (HRI has FWHM=0.2” and HST iLband-like fiuxes, LRI has FWHM=1.66” 
and IRACchl-like fiuxes), using three different methods: SExtractor dual-mode aperture, SExtractor dual-mode 
“best”, and t-phot. Vertical lines show the 5cr (dashed) and la (dotted) limits of the simulated LRI. Bottom: magni¬ 
fication of the top panel, showing only T-PHOT results, colour-coded as a function of the covariance index. See text for 
more details. 


4.1.4. Testing different fitting options: threshold fitting 

As described in Sect. |2.1.5| t-phot includes the option of 
imposing a lower threshold on the normalized fiuxes of tem¬ 
plates so as to exclude low signal-to-noise pixels from the 
fit. Figure p!5] shows a comparison of the relative errors ob¬ 
tained with three different values of the THRESHOLD param¬ 
eter: t = 0, t = 0.5, and t = 0.9 (whic means that only pix¬ 
els with normalized fiux fnorm > tx fpeak in the convolved 


template will be used in the fitting procedure). The differ¬ 
ences are quite small; however, a non-negligible global effect 
can be noticed: all sources tend to slightly decrease their 
measurement of fiux when using a threshold limit. This 
brings faint sources (generally overestimated without using 
the threshold) closer to their true value, at the same time 
making bright sources too faint. This effect deserves careful 
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investigation, which is beyond the scope of this study, and 
is postponed to a future paper. 


4.1.5. Colours 

A final test was run introducing realistic colours, i.e. as¬ 
signing fluxes to the sources in the LRI consiste nt wit h 
a realistic SED (as output by GenCat, see Sect. 4.1.2), 
instead of imposing them to be equal to the HRI fluxes. 
We took IRAC-chl as a reference filter for the LRI, con¬ 
sistently with the chosen FWHM of 1.66”. Furthermore, 
we allowed for variations in the bulge-to-disc ratios of the 
sources to take possible effects of colour gradients into ac¬ 
count. We compared the results obtained with t-phot with 
the ones obtained with two alternative methods to deter¬ 
mine the magnitudes of the sources in the LRI: namely, 
SExtractor dual mode aperture and MAG_BEST photome¬ 
try (with HRI as detection image). The differences between 
measured and input magnitudes in the LRI, iRmeas—^true^ 
are plotted in Fig. ^ Clearly, t-phot ensures the best 
results, with much less scatter in the measurements than 
both of the other two methods, and very few outliers. 


4.2. Performance on real datasets 


It is instructive to check how t-phot performs on real 
datasets, in addition to simulations. To this aim, we run 
two different tests. In the first, we compared the results 
of the TFIT CANDELS analy sis on the UDS CANDELS 
/-band (Galametz et al. 2013) to a t-phot run obtained 
using the cells-on-objects method and different parameters 
in the kernel registration stage. Figure shows the his¬ 
tograms of the differences in the photometric measurements 
between tfit and t-phot. Many sources end up with a 
substantially different flux, because of the two cited factors 
(a better kernel registration and the different fitting proce¬ 
dure). We note that the majority of the sources have fainter 
fluxes with respect to the previous meas ureme nts, precisely 
because of the effect described in Sect. 4.1.2[ fitting using 
a grid of cells introduces systematic errors assigning light 
from sources that are not listed in a given cell, but overlap 
with it to the objects recognized as belonging to the cell. 
To further check this point. Fig.shows some examples of 
the difference between the residuals obtained with tfit (of¬ 
ficial catalogue) and those obtained with this t-phot run 
using cells-on-objects method, also introducing better reg¬ 
istration parameters in the dance stage. Clearly, the results 
are substantially different, and many black spots (sources 
with spurious overestimated fluxes) have disappeared. Also, 
the registrations appear to be generally improved. 

The second test was run on FIR/sub-mm SCUBA- 
2 (450 /im, FWHM=7.5”) and Herschel (500 /im, 
FWHM=36”) images of the COSMOS-CANDELS field. In 
both cases, a list of 24+850 /im sources was used as un¬ 
resolved priors. Figure [TT] shows the original images in the 
top row, and the residuals in the bottom row. The model 
has removed all significant sources from the 450 /im map 
and the majority from the 500 /im map. Figure shows a 
comparison of the fluxes measured in the t-phot fits to the 
450 /im and 500 /im maps at 24+850 /im prior positions, 
with the error bars combining the errors on both flux mea¬ 
surements. Agreement within the errors implies successful 
deconfusion of the Herschel image to reproduce the fluxes 


measured in the higher resolution SCUBA-2 image. This 
typology of analysis is very complex and we do not want 
to addres s here the subtletie s of the process; we r efer the 
reader to Wang et al.| ( |2Q15[ in preparation) and Bourne 


eEliI] (|2Q15| in preparation) for detailed discussions on 
the definition of a robust and reliable approach to mea¬ 
sure FIR and sub-mm fluxes. These simple tests, however, 
clearly show that T-PHOT is successful at recovering the 
fluxes of target sources even in cases of extreme confusion 
and blending, within the accuracy limits of the method. 



Figure 17. Results from a test run using t-phot with 
unresolved priors on FIR/sub-mm real dataset. Upper 
row, left to right: SCUBA-2 450 /im (FWHM=7.5”) and 
Herschel 500 /im (FWHM=36”) images of the COSMOS- 
CANDELS fields. Lower row, left to right: residuals for the 
two fields, obtained with t-phot runs using a list of 24+850 
/im priors. See text for details. 


5. Computational times 

As anticipated, t-phot ensures a large saving of com¬ 
putational time compared to similar codes like tfit and 
CONVPHOT when used with identical input parameters. 
For example, a complete, double-pass run on the whole 
CANDELS UDS field at once (I band; ^35000 prior sources; 
LRI 30720x12800 ^ 400 million pixels; standard tfit pa¬ 
rameters and grid fitting) is completed without memory 
swaps in about 2 hours (i.e. 1 hour per pass) on a standard 
workstation (Intel i5, 3.20 GHz, RAM 8 Gb). A com- 
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UDS I band 




Figure 18. UDS I band tfit versus t-phot comparison. Top panel: compared measured fluxes. Bottom panel: histogram 
of relative measured flux difference. 



Figure 19. UDS I band tfit versus t-phot comparison. The panels on the left show two small patches of the official 
CANDELS residual image obtained using t-fit. The residual images of the same regions are showed in the right panels, 
this time obtained using t-phot with cells-on-objects method and improved local kernel registration. We note the 
disappearence of many spurious black spots. 


plete, double-pass run on the GOODS-S Hawk-I W1 field 
(^17500 prior sources, LRI 10700x10600 100 millions 

pixels, identical parameters) is completed in ^20 minutes. 
For comparison, TFIT may require many hours (^24) to 
complete a single pass on this Hawk-I field on the same 
machine. It must be said that tfit by default produces 
cutouts and templates for all the sources in the HRI image; 


selecting the ones belonging to the LRI field and inputting 
an ad hoc catalogue would have reduced the computing 
time by a factor of two (i.e. 11 hours for a single pass). 
It was not possible to process large images like the UDS 
field in a single run, because of RAM memory failure. CON- 
VPHOT timings and memory problems are similar to those 
of TFIT, although they have different causes (being writ- 
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measured Herschel 500 jam 


Figure 20. Accuracy of the flux determiuatiou for the dataset described iu Fig. 

(FWHM=36”) fluxes /500 are compared to the fluxes obtaiued for the SCUBA -2 450 /im (FWHM=7.5”) / 450 , cou- 
sidered as refereuce fluxes (usiug 24+850 /im uuresolved priors iu both cases). The symbols have the same meauiug as 
iu Fig. the error bars uow iuclude the measured error ou the refereuce flux. See text for more details. 


teu iu C, computatiou is geuerally faster, but it employs a 
slower couvolutiou method aud the solutiou of the liuear 
system iu performed as a siugle fit iustead of grid flttiug 
like iu TFIT, beiug much more time cous umiug ). 

Adoptiug the cells-on-objects (Sect. 2.1.5|) method iu- 
creases the computatioual time with respect to the tfit 
staudard cell approach, but it is still far more couveuieut 
thau the CONVPHOT staudard siugle-flt approach, aud gives 
uearly ideutical results. 

Table summarizes the computatioual times for ex- 
teuded tests ou a set of simulated images haviug differeut 
detectiou depths (aud therefore uumber of sources) aud di- 
meusious, with LRI FWHM=1.66”. The simulatious were 
ruu ou the same machiue described above, usiug three dif¬ 
fereut methods: whole image flttiug, cells-on-objects^ aud 
100 X 100 pixels cells flttiug. 


6. Summary and conclusions 

We have preseuted t-phot, a uew software package de¬ 
veloped withiu the astrodeep project, t-phot is a ro¬ 
bust aud versatile tool, aimed at the photometric aual- 
ysis of deep extragalactic fields at differeut waveleugths 
aud spatial resolutiou, decoufusiug bleuded sources iu low- 
resolutiou images. 

T-PHOT uses priors obtaiued from a high-resolutiou de¬ 
tectiou image to obtaiu uormalized templates at the lower 
resolutiou of a measuremeut image, aud miuimizes a 
problem to retrieve the multiplicative factor relative to each 
source, which is the searched quautity, i.e. the flux iu the 
LRI. The priors cau be either real cutouts from the HRI, or 
a list of positious to be fitted as PSF-shaped sources, or au- 
alytical 2 -D models, or a mix of the three types. Differeut 


optious for the flttiug stage are giveu, iucludiug a cells- 
on-objects method, which is computatioually efhcieut while 
yieldiug accurate results for relatively small FWHMs. T- 
PHOT eusures a large saviug of computatioual time as well 
as iucreased robustuess with respect to similar public codes 
like its direct predecessors tfit aud CONVPHOT. With au 
appropriate choice of the parameter settings, greater accu¬ 
racy is also achieved. 

As a fiual remark, it should be pointed out that the anal¬ 
ysis presented in this work deals with idealized situations, 
namely simulations or comparisons with the performances 
of other codes on real datasets. There are a number of sub¬ 
tle issues regarding complex aspects of the PSF-matching 
techinque, which become of crucial importance when work¬ 
ing on real data. A simple foretaste of such complexity can 
be obtained by considering the problem described in Sect. 

i.e. the correct amplitude to be assigned to the segmented 
area of a source. Work on this is ongoing, and the full dis¬ 
cussion will be presented in a subsequent companion paper. 

As we have shown, t-phot is an efficient tool for the 
photometric measurements of images on a very broad range 
of wavelengths, from UV to sub-mm, and is currently be¬ 
ing routinely used by the Astrodeep community to anal¬ 
yse data from different surveys (e.g. CANDELS, Frontier 
Fields, AEGIS). Its main advantages with respect to sim¬ 
ilar codes like tfit or CONVPHOT can be summarized as 
follows: 


— when used with the same parameter settings of tfit, 
T-PHOT is many times faster (up to hundreds of times), 
and the same can be said with respect to other similar 
codes (e.g., CONVPHOt); 
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Size[pix] 

27 

28 

29 

Number of sources 

2500 X 2500 

523 

1070 

1398 

5000 X 5000 

2104 

4237 

5561 

10000 X 10000 

8390 

16807 

22394 

20000 X 20000 

33853 

65536 

65536 

Whole image fitting 

2500 X 2500 

38” (2”) 

54” (10”) 

r9” (20”) 

5000 X 5000 

3’26” (I’l”) 

11’9” (7’41”) 

20’28” (16’1”) 

10000 X 10000 

lh28’22” (lhl5’46”) 

8h26’l” (71i58’10”) 

21hl6’24” (20h27’53”) 

20000 X 20000 

- 

- 

- 

Cells-on-objects fitting 

2500 X 2500 

46” (4”) 

I’ll” (16”) 

rso” (33”) 

5000 X 5000 

3’1” (18”) 

4’27” (1’8”) 

6’3” (2’20”) 

10000 X 10000 

12’27” (1’12”) 

17’52” (4’31”) 

25’11” (9’52”) 

20000 X 20000 

51’12” (6’1”) 

lh34’40” (35’8”) 

lh43’10” (41’2”) 

100x100 pixel 

s cells fitting 

2500 X 2500 

52” (3”) 

1’6” (7”) 

1T4” (9”) 

5000 X 5000 

3’16” (14”) 

4’22” (29”) 

4’54” (41”) 

10000 X 10000 

i3’4” be”) 

17’12” (1’54”) 

19’47” (2’53”) 

20000 X 20000 

55’24” (6’19”) 

lhl8’38” (15’53”) 

lhl7’17” (17’58”) 


Table 2. Test of computational times for t-phot runs on images of given dimensions and limiting magnitude in detection. 
Each entry of the table gives the total duration of run, the duration the fitting stage alone (in parentheses), and the 
number of fitted sources. The dance stage takes most of the CPU time after the fitting routine. 


— T-PHOT is more robust, more user-friendly, and can han¬ 
dle larger datasets thanks to an appropriate usage of the 
RAM; 

— T-PHOT can be used with three different types of priors 
(real high-resolution cutouts, analytical models and/or 
unresolved point sources) making it a versatile tool for 
the analysis of different datasets over a wide range of 
wavelengths from UV to sub-mm; 

— T-PHOT offers many options for performing the fit in 
different ways, and with an appropriate choice of pa¬ 
rameter settings it can give more accurate results. 

Future applications might include the processing of 
EUCLID and CCAT data. New releases of the software 
package, including further improvements and additional op¬ 
tions, are planned for the near future. 
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Appendix A: The parameter file 

Below is a template of the standard first-pass parameter file 
to be given as input to t-phot (similar templates for both 
the first and the second pass are included in the dowload- 
able tarball). It is very similar to the original tfit param¬ 
eter file, and part of the description is directly inherited 
from it. 

A.l. Pipeline 

Standard optical/NIR double-pass runs can be achieved by 
setting order standard and order standard2. 

A standard first-pass run includes the stages priors, 
convolve, fit, diags, dance, plotdance. The stage 
priors allows for an automatic re-construction of the 
pipeline depending on the input data given in the following 
sections (see the documentation included in the tarball). 
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# T-PHOT PARAMETER FILE 

# Control parameters: 



fitting coo 


# PIPELINE 


cellmask 

true 



mask!loor 

le-9 

# 1st pass 
order standard 


writecovar 

true 

#priors, convolve, fit, diags, dance, plotdance 

threshold 0.0 




linsyssolver lu 


# PRIORS STAGE 


clip false 


# Choose priors 

types in use: 

# DIAGNOSTICS STAGES 

usereal True 

usemodels True 


model!ile 

lores_collage_passl.fits 

useunresolved True 





# Dance: 


# Real 2-d profiles 

dzonesize 

100 

hires!ile 

HRI.fits 

maxshift 

1.0 

hirescat 

HRI.cat 

ddiagfile 

ddiags.txt 

hiresseg 

HRI.seg.fits 

dlogfile 

dlog.txt 

normalize 

true 

dancefft false 


subbckg True 
savecut true 

cutoutdir 

cutouts 



cutoutcat 

cutouts/_cutouts.cat 

A.2. Priors 



# Analytical 2-d models 
modelscat models/models.cat 
modelsdir models 


culling false 

# Unresolved point-like sources 
poscat pos.cat 

psffile psf.fits 

# CONVOLUTION STAGE 


loresfile 

loreserr 

errtype 

rmsconstant 

relscale 


LRI.fits 
LRI.rms.fits 
rms 
1 
1 


FFTconv true 
multikernels false 
kernel!ile kernel.fits 

kernellookup chl_dancecard.txt 


templatedir templates 
templatecat templates/_templates.cat 

# FITTING STAGE 


# Filenames: 
fitpars 
tphotcat 
tphotcell 
tphotcovar 


tpipe_tphot.param 
lores_tphot.cat_passl 
lores_tphot.cell_passl 
lores_tphot.covar_passl 


A standard second-pass run includes the stages convolve, 
fit, diags, archive. The archive stage creates a direc¬ 
tory after the name of the LRI, with some specifications, 
and archives the products of both runs. 

Double-pass runs for FIR/sub-mm can be achieved 
by setting order positions, fit, diags, dance, 
plotdance and order positions, fit, diags, 
archive. 


Each prior must have a unique identification number (ID) 
to avoid errors. The user must be careful to give the cor¬ 
rect information in this paramfile. Select the priors to be 
used by switching on/off the relative keywords: usereal, 
usemodels, useunresolved. 

— hires!ile: the high-resolution, detection image. If a 
catalogue and a segmentation map are given in the two 
subsequent entries (hirescat and hiresseg), cutouts 
will be created out of this image. This step is necessary 
if a catalogue of real or model priors are to be used. The 
catalogue hirescat must be in a standard format: id 
X y xmin ymin xmax ymax background SEx_flux (x 
and y are the coordinates of the source in HRI pixel 
reference frame; xmin, ymin, xmax, ymax are the lim¬ 
its of the segmentation relative to the source in HRI 
pixel reference frame; background is the value of the 
local background; and SEx_f lux is a reference isophotal 
flux). 

— poscat: a catalogue of positions for unresolved, point- 
like sources. No HRI image/segment at ion is needed, 
while the PSF to be used to create the models is manda¬ 
tory (psffile). The catalogue must be in the standard 
format id x y. 

— modelscat: a catalogue (with format id x y xmin 
ymin xmax ymax background SEx_f lux, as for a stan¬ 
dard HRI priors catalogue) of model priors, modelsdir 
is the directory in which the stamps of the models are 
stored. Models with two or more components can be 
processed, but each component must be treated as a 
separated object, with a different ID, and a catalogue 
for each component must be given. Catalogues for each 
component must have the same name, but ending with 
”_1”, ”_2”, etc.; put the ”_1” catalogue in the param¬ 
file. It is important to note that two components of the 
same object should not have exactly identical positions, 
to avoid numerical divergencies. 

— culling: if True, objects in the catalogue (real priors 
and/or models) but not falling into the LRI frame will 
not be processed; if it is false, all objects in the catalogue 
will be processed (useful for storing cutouts for future 
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reuse on different datasets) and the selection of objects 
will be done before the convolution stage. 

— subbckg: if True, subtract the value given in the input 
catalogue from each cutout stamp. 

— cutoutdir: the directory containing the cutouts. 

— cut out cat: the catalogue of the cutouts, containing the 
flux measured within the cutout area (which may be 
different from the SEx_flux given in the input cata¬ 
logue, e.g. if the segmentation has been dilated). We 
note that these are output parameters if you start from 
the priors/cutout stage; they are input parameters for 
the convolve stage. 

— normalize: determines whether the cutouts will be nor¬ 
malized or not; it is normally set to true, so that the final 
output catalogue will contain fluxes rather than colours. 

A.3. Convolution 

— loresfile, loreserr: the LRI and RMS images. T- 
PHOT is designed to work with an RMS map as the 
error map, but it will also accept a weight map, or a 
single constant value of the RMS from which an RMS 
map is generated. The errtype specifies which kind of 
error image is provided. For best results, use a source- 
weighted RMS map, to prevent the bright objects from 
dominating the fit. 

— relscale: the relative pixel scale between the two im¬ 
ages. For example if the HRI has a pixel scale of 0.1 
arcsec/pixel and the LRI has a pixel scale of 0.5 arc- 
sec/pixel, the value of relscale should be 5. If the LRI 
has been manipulated to match the HRI pixel scale and 
WCS data (e.g. using codes like SwARP by E. Bertin), 
put relscale 1. 

— kernelfile: the convolution kernel file. The kernel 
must be a FITS image on the same pixel scale as the 
high-resolutuion image. It should contain a centred, nor¬ 
malized image. 

— FFTconv is True if the convolution of cutouts with the 
smoothing kernel is to be done in Fourier space (via 
FFTW3). 

— kerntxt may be explicitely put True if one wishes to 
use a text file containing the kernel instead of a .fits 
one. T-PHOT supports the use of multiple kernels to ac¬ 
commodate a spatially varying PSF. To use this op¬ 
tion, set the multikernels value to true, and provide a 
kernellookup file (it is automatically produced during 
the dance stage in the first pass, but it can also be fed 
externally) that divides the LRI into rectangular zones, 
specified as pixei ranges, and provides a iocal convoiu- 
tion kernei fiiename for each zone. Any object that fails 
in a zone not included in the lookup file will use the 
transfer kernel specified as kernelfiie. 

— templatedir: the directory containing the tempiates 
created in the convoive stage, iisted in the cataiogue 
templatecat. We note that these are output parame¬ 
ters for the convoive stage, and an input parameter for 
aii subsequent stages. 

A.4. Fitting stage 

— fitpars, tphotcat, tphotcell, tphotcovar: these 
are aii output parameters. The tfitpars hie specifies 
the name of the speciai parameter hie for the fitting 


stage that wili be generated from the parameters in this 
file. The others are filenames for the output catalogue, 
cell, and covariance files, respectively. 

— fitting: this keyword tells T-PHOT which method to 
use to perform the fitting (see aiso Appendix B): 

— coo or 0 for cells-on-objects; 

— single or -1 for singie fit; 

— single! or -10 for optimized singie fit (the LRI 
is divided in square ceiis containing roughiy 10000 
sources each); 

— cell_xdim, cell_ydim, cell_overlap for an arbi¬ 
trary grid of ceiis. 

— cellmask: if true, it uses a mask to exciude pixels from 
the fit that do not contain a vaiue of at ieast maskf loor 
in at ieast one tempiate. 

— writecovar: if true, it writes the covariance informa¬ 
tion out to the tphotcovar hie. 

— threshold: forces the use of a threshold on the flux, so 
that only the central parts of the objects are used in the 
fitting process. 

— linsyssolver: the chosen solution method, i.e. LU, 
Cholesky, or Iterative Biconjugate Gradient (IBG). LU 
is default. 

— clip: tells whether to loop on the sources excluding 
negative solutions. 

A.5. Diagnostic stages 

— modelfile: the .fits file that will contain the collage 
made by multiplying each template by its best flux and 
dropping it into the right place. An additional diag¬ 
nostic file will be created: it will contain the difference 
image (LRI - modelfile). Its filename will be created 
by prepending resid_ to the modelfile. 

— dzonesize specifies the size of the rectangular zones 
over which the pixels’ cross-correlation between LRI and 
modelfile will be calculated during the dance stage. It 
should be comparable to the size over which misregis¬ 
tration should be roughly constant, but it must be large 
enough to contain enough objects to provide a good sig¬ 
nal to the cross-correlation. 

— maxshift specifies the maximum size of the x,y shift, 
in LRI pixel frame, that is considered to be valid. Any 
shift larger than this is considered spurious and dropped 
from the final results, and replaced by an interpolated 
value from the surrounding zones. Ideally, maxshift 
Ipixel X FWRMlri/FWEMhri. 

— ddiagf ile is an output parameter for the dance stage, 
and an input parameter for the plotdance stage. 

— dlogf ile is an output parameter; it simply contains the 
output from the cross-correlation process. 

— danceFFT: if True cross-correlation is to be performed 
using FFT techniques rather than in real pixel space. 

Appendix B: The cells-on-objects algorithm 

Experiments on simulated images (see Sect.j^ clearly show 
that fitting small regions (cells) of the LRI, as done by de¬ 
fault in TFIT, may potentially lead to large errors. This is 
particularly true if the dimensions of the cells are chosen to 
be smaller than an ideal size, which changes from case to 
case, but which should always be greater than ^10 times 
the FWHM. However, it can be mathematically shown that 
the “arbitrary cells” method intrinsically causes the intro- 
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duction of errors in the fit, as soon as a source is excluded 
from the cell (e.g. because its centre is outside the cell), but 
contributes with some flux in some of its pixels. 

Consider a cell containing N sources. For simplicity, as¬ 
sume that each source i only overlaps with the two neigh¬ 
bours i — 1 and i + 1. Furthermore, assume that a (A^ + l)-th 
source is contaminating the N-th source, but is excluded 
from the cell for some reason, for example (as in tfit) be¬ 
cause the centroid of the source lies outside the cell. 

The linear system for this cell AF = B will consist of 
a matrix A with only the elements on the diagonal and 
those with a ±1 offset as non-zero elements (a symmetric 
band matrix), and the vector B will contain the products of 
templates of each source with the real flux in t he LR I (as a 
summation on all pixels), as described in Sect. 2.1.4| Given 
the above assumptions, this means that the N-th term of B 
will be higher than it should be (because it is contaminated 
by the external source). 

Using the Cramer rule for the solution of squared linear 
systems, the flux for the object i is given by 


fi = 


det Ai 
det A 


(B.l) 


with Ai a square matrix in which the i-th columns is 
substituted with the vector B. If for example = 3, for 
i = 1 this gives 


fi = [Bi{A 22 Ass — ^ 23 ) “ ^ 12 (^ 2^23 — ^ 3 ^ 22 )]/detM 

(B.2) 

and since B^ is larger than it should be, fi will be overes¬ 
timated (slightly, if Ai 2 is not large, i.e. if sources 1 and 2 
do not strongly overlap). On the other hand, for i = 3 we 
have 


fs — [^ 11 (^ 22^3 — A 23 B 2 ) — ^ 12^3 + ^ 1 (^ 12 ^ 23 )]/detM 

(B.3) 

and in this case again might be small, but the first 
term given by ^ 11^22 will certainly be large, resulting in 
a catastrophic overestimation of /s. The value of /2 will of 
course be underestimated, as it would be easy to show. 

From this simple test case it is clear that arbitrarily di¬ 
viding the LRI into regions will always introduce errors 
(potentially non-negligible) in the fitting procedure, un¬ 
less some method for removing dangerous contaminating 
sources is devised. 

The cells-on-objects algorithm aims at ensuring the ac¬ 
curacy of the flux estimate while at the same time drasti¬ 
cally decreasing computationa l time s and memory require¬ 
ments. As explained in Sect. |2.1.4[ when this method is 
adopted a cell is centred around each detected source, and 
enlarged to include all its “potential” contaminant neigh¬ 
bours, and the contaminant of the contaminants, and so 
on. To avoid an infinite loop, the process of inclusion is 
interrupted when one of the following criteria is satisfied: 

— the flux of the new neighbour is lower than a given frac¬ 
tion ffiux of the flux of the central object (the consid¬ 
ered fluxes are: if real priors are used, the ones given 
in the HRI catalogue; if unresolved priors are used, the 
ones read in the pixels of the LRI containing the co¬ 
ordinates of the sources; if analytical models are used, 
the ones of the models as reported in the HRI models 
catalogue), or 


— the template of the neighbour overlaps with its direct 
previous contaminant for a fractional area lower than 

farea • 

Experiments on simulations have shown that good results 
are obtained with ffiux = 0-9 and farea = 0.25, and these 
values are used as constants in the source code. 

We note that if a cell is enlarged to more than 75% 
of the dimensions of the total LRI, t-phot automatically 
switches to the single fit on the whole image. 

Appendix C: Suggested best options 

Of course, different problems require different approaches 
in order to obtain their best possible solution, and users are 
encouraged to try different options and settings. However, 
some indicative guidelines for optimizing a run with T- 
PHOT can be summarized as follows. 

— Be sure that all the required input files exist and have 
correct format, and that paths are correctly given in the 
parameter file. 

— Whenever possible, fit the whole image at once (i.e. 
put fitting single in the parameter file). The more 
sources there are, and the more severe their blending, 
the more CPU time will be required (see Sect.|^. If the 
blending is not dramatic, it is safe to switch to the cells- 
on-objects method (i.e. put fitting coo in the param¬ 
eter file). On the other hand, if blending is severe this 
option would result in redundant fittings because cells 
would be enlarged to include as many neighbours as pos¬ 
sible, increasing the total computing time. In this case, 
either stick to the whole image fitting, or (depending on 
the desired degree of accuracy) switch to the TFlT-like 
cells fitting. 

— Spend some time in checking the output catalogue, e.g. 
considering with caution fits relative to sources having 
flags > 0 and covariance indices larger than 1. 
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